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THE ASYMPTOTIC EXPANSION OF THE STURM-LIOUVILLE 

FUNCTIONS. 

Bt F. H. Mtjerat. 

1. In the Proceedings of the National Academy of Sciences (vol. 3, 
pp. 656-659), Professor BirkhofT gave a direct proof of the closure of the 
set of Sturm-Liouville functions defined by the equation and boundary 
conditions, 

(1) S+^P'- ^(^)52/ = 0, 2/(0) = yd) = 0. 

This proof was based on a certain asymptotic expansion for these 
functions, similar to those used by Hobson* and Kneser.f At the sug- 
gestion of Professor Birkhoff I have undertaken to develop this expansion 
in detail, using explicitly the Volterra integral equation of the second kind 
used more or less implicitly by most writers in this connection; the method 
of successive approximation employed in the asymptotic development of 
the characteristic numbers is capable of extension to the functions satis- 
fying the boundary conditions 

a'tjiO) - at/'(0) = 0, aa' ^ 0, 

^'j/(0) + ^2/'(0) = 0, 00' ^. 0. 

The explicit use of the Volterra integral equation is especially con- 
venient in the study of the differentiability of the characteristic functions 
with respect to a parameter a, introduced by replacing g{x) by agix). 

1. Some preliminary inequalities. Assume that g(x) has bounded vari- 
ation; instead of the system (1) consider first the system 

(2) S + [''' - <^^^^)5^ = 0. 2/(0) = 2/(1) = 0. 
The equation above can be written in the form 

TO 

^ -f p-y = <jg{x)y, 
which leads to the Volterra integral equation of the second kind, 

* Proceedings of the London Math. .Soc, ser. (2), vol. 6, p. 374. 

t Die Integralgleiehungcn und ihre Anwendung in der mathematischen Physik, Chap. 3. 
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(3) y(x) = a sin p{x - /3) + ^ fsin p(x - Og(.Oyi^)d^- 

P Jo 

Here a. 8 are arbitrary constants; if y(x) satisfies the boundary condition 
j/(0) = 0, we may assume /3 = 0. For convenience assume a positive 
or zero; p, ff are real, and g(x) is real for :S x :^ 1. Since g{x) has 
bounded variation ior :^ x ^ 1,\ g(x)\ has an upper bound G in this 
interval. 
Suppose 

Z.Ca-il) = sin p(x - Ogi^), I K^{x,0 I ^ <?• 
Equation (3) becomes, 

(4) y(x) = a sin px+-f K^(z, Oy(^)d^. 

P Jo 

Suppose 



^^'=(^yrr---r"'^^^^'^^^^^^^'-'^) 



(5) 

• • • -K,(^,_i, Q sin pMlid^2 • • • <i^.- 

Then the solution of (4) can be given in the form 
(6) y(x) = a I sin px + Z l\{x) \ , 

and this series is dominated by the series 

For convenience assume a = 1. From (5), (6) it is seen immediately 
that for large values of j p | the first few terms of the series are the impor- 
tant ones. Expand Uii 

Ui = - I sin pix - ^i) sin p|i?(^i)d$i. 

P Jo 

Substitute 

sin p^i sin p(x — $i) = §[cos p(x — 2|i) — cos px], 

Ui = -^cos pxj g{ki)d^i +YpJ ^^^ ^^^ ~ -^i)S'*^^i)'^^i- 

From a lemma by Riemann* it follows that the second term on the 
right is of the order of 1/p-. 



' Gesammelte Werke, p. 241; Whittaker and Watson, Modern Analysis, 2d ed., p. 166. 
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(8) 



Jo 

H{x, a, p) = fp- J'cos p(x - 2i,)g{^,)di, + gC^,. 
Then from (6), 



(9) 



y{x) = sin px — 



(Th(x) 



cos px + H(X, <T, p). 



The function y(x) has already been so determined as to satisfy the 
first boundary condition t/(0) = 0, for all values of p; it remains to deter- 
mine the particular values of p for which y{l) = 0, or the characteristic 
numbers. This condition becomes: 



(10) 



tan p = -^ = <t>((r, p). 

2p cos p -^v ' f/ 



This equation can be solved by a method of successive approximations 
involving a Lipschitz condition of the form 

I </.(<r, p") -<A(<r, p')| <C,\p" - p'\. 

To find an upper bound for C,, compute the partial derivatives dUJdp: 



dU, 
dp 



Since 



= -lU.ix,.,p)+{iy££...£-[tK^ix,,.) 

■ ■ ■ K,{^k-„ a-i, P) ^ i^p(l*-i, ?*)••• -Kp(^.-i, Q ] [sin pU^id^2 ■ ■ ■ d^,] 



^A%(|,_„ IJ = ($,_, - tj cos p(^,_i - $J6r(?,), 



Consequently 
d 



^x,a._:, ui^G. 



(11) 



C7,(x, <r, p) 



(»<■ 



(12) 



ap 

From (8), 

air (T C 

= -nT— cos p(x - 2ti)sf(|i)d^i 

-P' Jo 



:h),(fnf('+->^)] 



dp 



- :f r (x - 2|i) sin p(x - 2,'i)(7(ti)d?i + Z 

^P Jo ^=:: 



dp 
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It follows from (11), (12), and Riemann's lemma that there exists a 
constant Oi such that if p is real and greater than 1, 



(13) 



From (10), 



(14) 



dp '^p-' 

d4> -<rh eosp^ + Fsmp 



dp 2p- cos^ p 

From (7), (8), a constant Oo can be found such that if p > 1, 



(15) 



H(X, (7, p) I < -^, 



If p is so chosen that ] cos p | s §, and \ p\ > 1, it follows from (13), 
(14), (15) that for some m, 



(16) 



dp 



m 
p- 



This inequality holds forO:^x^l,0:S<r^l, and p real; it leads 
immediateh' to the Lipschitz condition desired. 

It will be convenient also to calculate d4>;d<j; from (10), 

dH 

d4> h dcr 

d(T 2p cos p * 



If I cos p I > 5 

U<r~ 2p 



< 2 



\dH\ 
do- 



Since from (5), 

\ dU, 
I dcr 

we have finally, 

(17) 



< 2^ fcos p(x - 2t0g($i)d,'i + Z^ 
- p'\p) A - p\ p ) ' {v-\)\' 



\ d4) h_\ ^ Oi 
\d'(r~2'p\ ~7' 

where as is a constant; p is real and \p\ > 1, | cos p | 2: 5. 

2. Asymptotic expansion of the characteristic numbers. In the equation 

tan p = </>(<r, p) 

put p — Kir + i,; the equation becomes 

(18) tan e, = 4>{cr, kt + «,). 
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Let 5, be the interval kt — (tt 3) ^ p ^ kit + (^,3). From the 
inequalities already established it is seen immediately that N can be 
chosen so large that if « > JV, and if pi, po are in the same interval 5,, the 
following inequalities are satisfied : 

c 

(a) I (/>((r, pi)| < , _ ^, , 

rh-\ \ ^f \ j.( \\ ^ w I P2 - Pi I 
(0) I <^(o-, P2) - <^(<r, pi) I < r/ _ Y^^Y ' 

(19) (c) ^^ < 1 

, ,, 2C 1 TT 

^''^ {k - l)7r ■ 4m ~ ^ 6 • 

[(k - l):rP 

Condition (b) follows from (16) by integration. Let e', «", • • ■, €<"' 
be defined by the equations 

tan «' = 4>((T, kt), 

(20) ^^^ *" = "^("^j "'^ + *'). 

tan €(■■> = <i>{(T, Kir + e^-"). 

In solving this system, observe that if | m i < 1; 

(21) arctan /x = I :; — , — ; , arctan u < r-^ — —^ ■ 

'^ Jo 1 + X- ' ' '^ ' 1 - M^ 

From this inequaUty we obtain : 

I , , I <t>{<^, Kir) I 



1 — ^(c, kt)- ' 

From (19); (a) and (d), 

I 4>(<T, Kir)\ <j^. 

Consequently 

I <i>-{(T, ktt) I < 5, 

(22) ! e'l < 2| 0((r, ktt) | <^. 

It follows that {kit + e') lies in the interval 5,, and e" can be determined 
from (20). 

Applying the method of induction, we shall assume that Uir + e'), 
{kit + i") • ■ •(kh- + €<">) He in the same interval 5.; in addition, assume 
that for M < «', 
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^ ' ' * ' ^ (/c - 1)^L[(« - 1)7^]= J ' 

tan [.<^+^) - e<^)l = <^(^, K^ + ^'-') j-,^^_r^, '--^ + «'-'jO_ 

From (19) (5), (c?), and from (23), 

I tan fe<^+» - e«ll < ^^ . ^2C_ I" __4m_ T- 
,tanLe 6 Ji<[(, _i)^]2 (, _i)^[[(, _l)^]2j • 

From (19) (c), (d), and from (21), 

r,(.^i) _ ,(.)i < — g£- r 4m -y 

^' ' J ^ (k-DttLkk-d^pJ' 

I fC+i) I = i e' + (e" -«')+•••+ (^^'^^ - e<"') i 

2C r 4m r 4w in 

(24) <(«- l)xL^^[(/c-l)^P+ ■■■ ^LK^- 1)t?J J 

TT 

from (19), (d). Consequently the inequalities (23), (24) hold for y = 1, 
and «'" can be determined from (20), while {k-k + «'), (^'t + e") lie in 6,; 
[kv + e'") lies in 5„ hence 6*''* can be determined, and (ktt + e'*^) lies 
in 5,, etc. By the principle of induction (23) and (24) hold for i- = 1, 2, 3, 
•••,!', • • • ; from (23) the series 

e. = 6' + (e" - O + («'" - e") + • • • 

converges absolutely and uniformly with respect to cr, for :£ <r :S 1; 
consequently e, is a continuous function of c 

To verify that e, satisfies (18), observe that in the identity 

tan €, — <^((7, KTT + e,) = [tan e, — tan e^'+^'j 

+ [tan fC+i' - (A(<r, KTT + f^'O] 

+ [<^((r, Kir + e(^^) - (j!.((r, ktt + €j], 

the second term on the right vanishes, while the first and last terms can 
be shown to approach zero, as v becomes infinite, with the aid of the 
inequalities (19) and (23). Since the difference on the left is independent 
of V, this difference must vanish. 

The asymptotic expansion of «, can be determined with the aid of the 
inequalities established above. 
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I e, - € 1 = |(e -«) + (« -«)+••• I 
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2C r 4m "I 



< 



47n 



6 [(k - 1)tP' 
e' = arctan (^(<r, ktt) 

= 0((r, ktt) + C3[()!>((r, KTr)]^ ± 
hu H 



e, — 



2K7r cos Kir 

2/C71 



+ 



<-<■)+ (<■-£)! 



(25) 



_ hff_ H((T, k) 
Zktt k-tt- 



where H is a function less than some constant H for all values of « > A^, 
and for :s 0- ^ 1. For large values of k the derivative dp, I da exists. 
For suppose ci, o-o two values of a in the interval ^ o- ^ 1, pi, p2 the 
corresponding values of p for a given k. 
From the equations 

tan pi = 4>{cT\, pi), tan pj = <i>{(j2, pz), 



tan P2 — tan pi r, , , , , tan (p2 — pi) 

= [1 + tan pi tan p2J 



we obtain 

P i — Pi . 
<T2 — <^l 



P2 — Pi 



P2 — Pl 



<A(<T2, Pl) - <t>(cri, Pl) 



(To — (Ti 



[1 + 4>icri, Pi)<t>{<T2, Pi)] 



tan (p2 — Pl) (f>((T2, P'i) — 4>(a-i, pi) " 



P2 — Pl P2 — Pl 

Let the difference a^ — <ri approach zero; the partial derivatives ap- 
proached on the right exist and are continuous in p, o-; consequently 



dp 
da 



30 
da 



op 



30 
da 



1 + 



30 
3p 



- <t>"- 



1 +.0= - 



30 

dp 
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Consequently from (11), 

(96) ^ = A + EM 

^-^^ da 2p^ p- ' 

where ivp is a function of p, a, bounded for large values of p. 
3. The characteristic functions. In equation (9), 

crhlx) 
(9) y(x) = sin px ^ — ^^s px + H(x, c, p), 

substitute p = ktt -\- e^: 

sin px = sin kttx + e.x cos kvx + e,-i^i(x, k, c), 

cos px = cos KirX + e,i/'2(^, K, f), 

where ^\, rp-i and hereafter ^p^ will denote functions of x, k, a continuous 
in X and less than some constant C„ for k > X, and :^ o- ^ 1. 
Equation (9) becomes, 

(9') y(x) = sin kttx + ?, — [xh — /iCar)] cos kttx + — . 

2 Kir K- 

Equations (9) and (9') are identical for ^ <r S 1. In -(9) com- 
pute dy/da: 

dy dH h(x) 

(20 +^^x cos p:r + -3^ cos px + -27-sm P^ + ^ J 

xh — h{x) , ;/'4 

= H cos KirX H 17 . 

/kit k- 

From (9'), 

9y xh — h(x) , 1 d\(/i 

-^ = — o cos KTTX + — -3— . 

Off ^Kir K- acr 

Comparing results, 

dxps 

From (9') compute 

I y-(x)da 

Jo 

I y-(x)dx = I sin- kxx + x — [xh — h{x)] sin 2Kirx -\ — ^ \dx 

(2S) 

1 , 'As 



XX 



= + 



.-2 ' 
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since from Riemann's Lemma, 



Also, 



I n — i^h — h(x)] sin 2Kirxdx < — . 
j^j^ y'-ix)dx = 2 J^ i/(a;) I 2^(a;/i - /i(x)) cos ktx + -^ \dx, 



applying Riemann's Lemma again, we obtain the result 



A 
do 



c' 
<T2> 



where c' as well as c is a constant independent of k if k > A'. From (28) 
(29) 



34^ 



< C , ^ = -At. 



Each of the functions j/,(x) is normalized by division by the function, 

of (T, 



"Vj 2/<^Ca;)d-r. 



= 2 - 



V^=V2 + V2[Vrrf.-l] 






If !2i < 1, 



.,—— , 3 1-1 , , 11-3 3^ 



Vl +2 - 1| < V[l +1^1 +|2l'+ •••] 



< 



2 

1 I 2 



2 1 - 1 2r 

If the constant A" has been chosen sufficiently large, 



and accordingly 



Hi 

\2k- 



1 

<0' 



J2+^!=V2+^. 

V K- K- 
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(30) Vi 2/'^^)^^ 

= V2 •[ sin (cttx + ^ — [x/i — h(x)] cos kttx H — t \ . 
[ Zktt k- \ 



(31) d<7 



j y-{x)dx 



r^\xh - h(x) l/'u 

= ^2 — 7s cos Kirx -\ 5- 

Zktt k- 



Comparing (30), (31), we obtain the result: 
(32) -j^ = ^n. 

Hence, finally, the normalized functions j/<(x) satisfying the sj-stem (2) 
are given asymptotically in the form 

y,(x) = V2 sin kitx + kt {^^'^ — ^(^)] cos kvx H —^ — \, 

bi/{x, <r, k) I 



(32) 



3(7 j 



<K, 



where K is independent of k, <r < 1 for k > iV. 

4. Nature of the solutions for small values of p. In the preceding discus- 
sion the parameter p was assumed larger than a constant AT'; to complete 
the discussion it is necessary to study the functions y{x) for small values 
of p. For this purpose write the equation of (2) in the form 

(20 % = \<jg{x) - \\y. 

The solutions of the homogeneous equation obtained by equating the 
left-hand member to zero are 1 and x; the corresponding integral equa- 
tion becomes, 

2/(x) = a + /3x + r(x - ^)\<jg{i) - \]y(^)d^, 
since y{0) = 0, a = 0; for convenience assume /3 = 1, 
y{x) = x+ Cix - ^)[<Tg{^) - \]y(M)dl 

Jo 
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For values of X less than or equal to N'' the kernel is bounded, and con- 
tinuous and different iable with respect to o-, X; the series expansion of the 
solution of this integral equation converges uniformly as before for 
:^ <r :£ 1, | X| ^ N''. Consequently the solution y{x) is continuous 
in X, (7, and as before the partial derivatives with respect to these par- 
ameters can be computed from the series expansion. 

If dy(d\ has continuous first and second partial derivatives with 
respect to x, such that 

dx\d\) d\\dx)' 



r- {dy\^±fd:y\ 
x-\d\J d\\dx^J' 



dx- 
then from the equation 

(2') y" + (X - a9{x))y = 0, 

we obtain: 

(33) S(l!) + (^-''W)S = -»- 

Let y(x) be a solution of (2') vanishing for x = 0, and suppose ^(x) 
any solution not vanishing at the origin for any value of a-. The solution 
of (33) becomes, 



dy 



= cjy(x) + c,y{x) + ^-, 1 -y, £ [y(x)y{^) - y{x)m][- vi^d^- 
When X = 0, dy!d\ = 0, y{0) = 0, hence Ct = 0. 



dy 



= cyix) - , ^ , r[ij{x)yiO - yix)ymy{^)d^. 
yy yy Jo 



The conditions d'-y dxd\ = d'-y d\dx, etc., are easily seen to be satisfied. 
Suppose X a characteristic number; then j/(l) = 0. Put x = I. 

[dyl - yd) f' ,,,,.,, 

Since 2/(1) = 0, y{l) 4= 0, hence 

+ 0. 

It follows immediately that dX/dv can be computed from the equation 

Vdyl rdyd_Xl ^ 

and is finite for p < N'. 
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Consequently if <ti, cto lie in the interval s a- < 1, cri < crj, 

< --i, ! 0-2 - 0-1 I, 

where A, is independent of cr, for :^ o- < 1. 

Since N' is finite, some A, is equal to or greater than any of the others; 
call this one A : 

(34) j y,{x, a.) - yjx, cri)\ < A\<T2 - ffi\ 

P. < N', S CTi < (72 S 1. 

It has been seen that y^ix, cr) is continuous in <r; it follows that 2/,(x, cr) 
vanishes just (k — 1) times in the interval < a; < 1. For this is true 
when (T = 0; as (T increases from to 1 the end-points of the curve y^yjx) 
remain fixed, and consequently zeros can enter or disappear only if for at 
least one value of cr the curve becomes tangent to the j-axis. At the 
point of tangency y =^ y' — 0; since y{x) satisfies the differential equa- 
tion (2) yix) vanishes identically. But from (32) this is seen to be 
impossible for very large values of k; suppose N so large that for k ^ N 
yJx, a) does not vanish identically for any value of cr between and 1. 
Then y,{x, a) has the same number of zeros for cr = 1 as for a- = 0, or 
(*c — 1) zeros. Now for k ^ A' the characteristic functions can be so 
ordered that yj,x, 1) has one more zero than 2/,_i(a;, 1); assuming the 
functions j/, in this order, j/.v-i has (A" — 2) zeros, y^^t has {N — 3) 
zeros, — hence for k ^ N a,s well as for k > N, y, has just (k — 1) zeros 
in the interval 

< a; < 1. 

The equation and inequality (32), together with (34), are sufficient 
for an immediate application of the theorem of Professor Birkhoff, with 
the aid of which the closure of the set of normalized functions satisfying 
(1) is established. 



